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ABSTRACT 

We investigate the non-gravitational heating of hot gas in clusters of galaxies (in- 
tracluster medium; ICM) on the assumption that the gas is heated well before cluster 
formation ('preheating'). We examine the jet activities of radio galaxies as the sources 
of excess energy in ICM, and the deformation of the cosmic microwave background (the 
Sunyaev-Zel'dovich effect) by hot electrons produced at the jet terminal shocks. We 
show that the observed excess entropy of ICM and COBE/FIRAS upper limit for the 
Compton y-parameter are compatible with each other only when the heating by the jets 
occurred at relatively small redshift {z < 3). Since this result contradicts the assump- 
tion of 'preheating', it suggests that the heating occurred simultaneously with or after 
cluster formation. 

Subject headings: galaxies:clusters:general-intergalactic medium-galaxies: jets-cosmic 
microwave background-X-ray : galaxies : clusters 

1. Introduction 

The departure of the properties of X-ray emitting gas in galaxy clusters (intracluster medium; 
ICM) from simple scaling relations gives rise to arguments about its thermal history (Kaiser 1986, 
1991; Evrad & Henry 1991; Fujita & Takahara 2000). Observed relations between X-ray luminosity 
and temperature show that from a rich cluster scale to a poor cluster scale, the exponent increases 
from Lx oc r|"^ (e.g. David et al. 1993; Xue & Wu 2000) to Lx oc r| (Ponman et al. 1996; Xue & 
Wu 2000), which are much steeper than those obtained by gravitational collapse alone {Lx oc T^-). 
Moreover, the discovery of the excess entropy in poor clusters ("entropy-floor") by Ponman et al. 
(1999) is presently interpreted as a strong evidence for the existence of non-gravitational heating 
in the ICM. 
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One of the popular scenarios to successfully explain these thermal properties of ICM is a 

so-called "preheating" scenario (e.g. Tozzi et al. 1999). In this scenario, it is assumed that the 
proto-ICM is heated well before the collapse of clusters, which explains the break in the Lx — Tx 
relation and the entropy floor as follows. If the virial temperature of a cluster, Tyir, is much higher 
than the temperature of the external gas, Tgx, the external gas accreted by the cluster is heated to 
Tvir by shock waves forming at the collapse, and the scaling relation Lx oc should be satisfied 
(e.g. Cavalicrc, Mcnci & Tozzi 1998). On the other hand, if Tyir is comparable to or smaller than 
Tex, shock waves do not form at the collapse, and the external gas adiabatically accretes onto the 
cluster. Thus, the gas temperature is determined not only by T^i^ but also by Tex- Several authors 
have shown that preheating models can reproduce the observational results (e.g. Tozzi et al. 1999), 
however, the heat source itself and the input epoch have not been identified. 

Some authors have investigated the heating by supernovae. However, Valageas &; Silk (1999) 
showed that the energy provided by supernovae cannot raise the entropy of intergalactic medium 
(IGM) up to the level required by current observations. Moreover, Kravtsov & Yepes (2000) 
estimated the energy provided by supernovae from the observed metal abundance of ICM and 
found that the heating by supernovae alone requires unrcalistically high efficiency. On the other 
hand, AGNs may be much more powerful and are, therefore, plausible candidates of heating sources, 
and thus we focus on AGNs. 

As for the input epoch, there have been few concrete arguments except for the ones inherent in 
models of individual sources. For AGNs, there have been no additional constraints like the metal 
abundance in the case of supernova heating (Kaiser &; Alexander 1999). In this paper, we propose 
a new approach to study this subject. Preheating sources which are so powerful as to influence 
the thermal properties of ICM would also deform the spectra of the cosmic microwave background 
(CMB) via inverse-Compton scattering (Sunyaev-Zel'dovich effect; hereafter SZ). Indeed, it has 
been shown that the energy supplied by AGN jets could be a significant source of SZ effect (Yamada 
et al. 1999). If the preheating scenario is correct, then the ICM should have been heated before the 
collapse of poor clusters. This allows us to obtain a lower limit of the redshift of energy input. On 
the other hand, if a large amount of hot IGM exists too early, the cumulative SZ effect would break 
the constraint by COBE/FIRAS (y < 1.5 x 10"^, Fixen et al. 1996). In this paper, we compare the 
estimated SZ effect and the excess energy induced by cocoons formed by AGN jet activities with 
the observational constraints, and discuss the validity of the simple preheating scenario. 



2. Models 

2.1. Cocoon model and the Compton y-parameter 

We assume for simplicity that the heating of proto-ICM by jet activities occurred well before 
cluster formation. The kinetic energy of the jet is transferred to thermal energy of the proto-ICM 
via thermalization at the shock at the hot spot. The thermalized jet matter expands into the 
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intergalactic medium (or proto-ICM) surrounding the radio galaxy laterally as well as along the 
jet axis (Begelman & Cioffi 1989; Nath 1995; Kaiser k Alexander 1997; Yamada et al. 1999). As 
this hot matter expands supersonically, a hot region surrounded by a shock surface around the 
radio galaxy is expected to form; hereafter we refer to it as a "cocoon". We briefly summarize the 
evolution of the cocoon below (see for details, Yamada et al. 1999). 

While the jet is active, we can write pressure inside the cocoon as 



by balancing the shock thrust and ram-pressure of the background matter, and by using Rankine- 
Hugoniot conditions. In this equation Lj represents the kinetic energy of the jet, 7 = 5/3 is the 
adiabatic index, is the "volume factor" which describes the shape of the cocoon (compared with 
a sphere), t is the time elapsed since the ignition of the jet, (p is the opening angle of the bow shock 
ahead of the jet termination spot (see Fig. 1 of Yamada et al. 1999), Ca is the sound speed, pa is 
the gas density. Pa is the pressure, and the suffix a denotes the ambient matter, respectively. In 
the derivation of the above equation, we assumed that Lj is time-independent, and is given by the 
Eddington luminosity of the central black hole which activates the jet. We use the gas density of 
the background universe as pn- ^Ve adopt the black hole mass -Mbh — 0.002A/gpjj, where Mgp}^ is 
the mass of the spheroidal component of the host galaxy (Magorrian et al. 1998). We can estimate 
the total thermal energy deposited during the jet active phase as PcV = eg x (7 — l)Ljtiife, where 
V is the volume of the cocoon, eg is the thermalization efficiency, and tiHe is the lifetime of the jet 
activity. We adopt the standard value for tute (3 x 10^ years), and take the thermalization efficiency 
for a free parameter. 

While the jet lifetime is short, the cocoon can stay hot for about its cooling time even after 
the energy supply stops. This "residual" phase contributes much more strongly to the SZ effect 
(Yamada et al. 1999). We model the evolution of the cocoon adopting the analogy with the evolution 
of a supernova remnant (SNR) in interstellar medium. Thus we define the effective lifetime of the 
cocoon as the time from the death of the jet to the epoch when the expansion time Tex = Rs/vg 
eqTials the cooling time behind the shock front, Tcooilshock {Rs Sind Vs are the shock radius and the 
velocity, respectively) . 

Finally we write the evolution of the internal energy of a single cocoon: 



where t is the cosmological time, tbegin is the epoch of jet ignition, and tr{z) is the effective lifetime 
at redshift z, respectively. Numerical simulations have shown that, although the thermal energy 
of a SNR rapidly decreases soon after the radiative phase begins, the cooling time increases as the 
SNR expands further and about 10% of the initial thermal energy is left behind (Chevalier 1974; 
Thornton et al. 1998). Thus we keep the total internal energy constant after it drops to 10% of its 
initial value. 




(1) 




(2) 
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The Compton y-parameter is calculated by integrating the product of total internal energy 
within a single cocoon and the number of radio galaxies in the line of sight, 

f f P V 1 

—^crTnRGiM,Zcoii)a^r'^dr—^dM, (3) 

where rUe is the electron mass, gt is the Thomson scattering cross section, riRQ is the comoving 
number density of radio galaxies, ZcoW is the typical collapse epoch of host galaxy halos, a is the 
cosmological scale factor, r is the comoving radial coordinate, Ra is the angular diameter distance, 
respectively (Yamada et al. 1999). We assume that radio galaxies reside in halos with the mass 
of M > 10^° Mq. We also assume that a fraction of normal galaxies has jet activity, and set the 
proportional constant to be a canonical value fr = 0.01. We count the number of radio (or normal) 
galaxies using the Press-Schechter number density (npg), and do not use the luminosity function 
of radio galaxies; this is because we intend to count up the "residual" cocoons, whose effective 
lifetime {tr) is much longer than the synchrotron-decay time. We define the epoch of jet ignition 
tbegin separately from the "typical collapse time" tcoU of the dark halo of the host galaxy, at which 
the variance of density perturbation of scale M, is equal to 1.69^ when f^o = 1. We take 

igap = ibegin ~ ^coU ^ & ^ee parameter because of the uncertainty inherent in the jet activation 
mechanism. The value of ^gap varies between (jet ignition coeval with the collapse of the dark 
halo of its host galaxy) and < 10^° years {k. Hq^). 



2.2. Energy Input into Proto-ICM 

It is reasonable to assume that the density of proto-ICM traces that of galaxies before the 
density perturbation corresponding to a protocluster (5) goes nonlinear. Thus the local number 
density of radio galaxies and the gas density of proto-ICM in a proto-cluster region are written as: 

%ai(M,2;) = nps/.(^)'(l + <5), (4) 

ngas(M,z) = P^n^f,{l + 5), (5) 
jinip 

where ag is the present scale factor, p^rit is the critical density, /i = 0.59 is the mean molecular 
weight for primordial gas, nip is the proton mass, and /c is the fraction of the gas compared with 
the baryon density f^b, respectively. 

According to equation (4), the energy density ejected by AGNs into a proto-cluster is given as 
etot = / /rnps(M, z) ■ PcV{M, z) ■ {1 + 5)dM. (6) 
Hence the energy input per nucleon E^^^^^ at present {z = 0) is 

£^input = = (^) ' 7^ / nps(M, 0)/, • P,V{M, 0)dM, (7) 
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which measures the additional, non-gravitational heating. Hereafter we assume /c = 1, which is 
reasonable when the density contrast is in a linear regime. 

3. Results 

We calculate equations (3) and (7) with various combinations of two parameters eg and 
tgap- We adopt the standard cold dark matter cosmology, with = l,h = 0.8, baryon den- 
sity Oft/i^ = 0.0125, and COBE normalization for the density perturbations. In Figure 1 contours 
of the Compton y-parameter and -Einput are plotted. X-ray observations show that excess energy 
via non-gravitational heating is 0.44ib 0.3 keV per a particle (Lloyd-Davies, Ponman &; Cannon 
2000). The region where parameter values are consistent with the y-parameter constraint obtained 
by COBE, y < 1.5 X 10"^ (Fixen et al. 1996), and the energy deduced by Lloyd-Davies et al. (2000) 
is indicated as a shaded region. As is clearly seen, almost horizontal contours from the y-parameter 
constraint severely limit the value of igap to be > 6.3 x 10® years. 

The result that the value of y is almost independent of thermalization efficiency eg comes 
from the weak dependence of the expansion speed on the internal energy in the "Sedov" phase 
{vs oc E^^^, Shu 1992), which results in the weak dependence of cooling time (effective lifetime of a 
cocoon) on eg. On the other hand, contours of -Einput limit the thermalization efficiency to eg < 0.4. 
The fact that the contours of -Einput run vertically for small values of tg^p reflects the rapid cooling 
of cocoons at high redshift. When tgap is small, the cooling time of a cocoon is so short that the 
cocoon energy rapidly reduces to 10% of the energy supplied by the jet; thus -Einput is simply given 
by Ljtiifeeff x 0.1 and is proportional to eg. For late input (large igap) the contours slightly curve to 
left; this is due to the small number of cocoons that contribute to the heating, which means that a 
large value of eg is needed to gain the same amount of energy. In order to assign tgap and the jet 
ignition epoch, we plot the corresponding redshift ^begin as the function of halo mass in Figure 2. 
Figure 1 shows that tgap ^ 6.3 x 10® years, which means heating of ICM occurred at Zbcgin ^ 3 
(Figure 2). Compared to the halo collapse epoch (tgap = 0; solid line), the upper bound is very 
close to the formation epoch of poor clusters with the mass of M = 10^3 _ iqI^Mq. This suggests 
that the heating of the ICM occurred simultaneously with or after the collapse of the poor clusters, 
which contradicts the assumption about the heating epoch inherent in the scenario. In other words, 
our results bring up a serious question about the genuine preheating model. 

4. Discussion and Conclusions 

We have proposed a new way to elucidate the thermal history of ICM. We have calculated the 
Compton y-parameter and the energy ejected into the ICM through AGN jet activities assuming 
that the non-gravitational heat input by the jets occurred well before cluster formation. Comparing 
them with observations, we have found that the heat input had not occurred at z > 3. Since 
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z ~ 3 is the typical formation epoch of poor clusters, this is not consistent with the assumption of 
'preheating' and suggests that the heat input occurred simultaneously with or after the formation of 
poor clusters. Considering the wide range dispersion in the formation epochs of poor clusters (z < 3; 
Lacey & Cole 1993; Kitayama &: Suto 1996; Balogh et al. 1999), the scenario of heating coeval with 
the cluster formation may be more plausible and would be consistent with the dispersions of cluster 
properties (Fujita &; Takahara 2000). Note that if the heating occurred in dense environment like 
in clusters, the lifetime and filling- factor of cocoons will decrease (Yamada et al. 1999), which 
may reduce the expected value of the y-parameter well below the observational limit. Valageas & 
Silk (1999) also estimated the y-parameter for QSO heating, but found a smaller value than ours 
{y ^ 10~^). This may be because in their model, the energy injected into IGM by QSOs started to 
cool immediately after the energy injection, which leads to a smaller y value. On the other hand, 
in our model, a cocoon remains hot for a long time until the cooling sets in (see §2). 

Below we discuss several points which are omitted in our simple model. First, a part of 
electrons may be accelerated at the terminal shock to become non-thermal populations. The SZ 

effect concerning these populations was calculated by several authors (see e.g., Birkinshaw 1999; 
Ensslin & Kaiser 2000), and is shown that the amplitude of the signal is reduced by only a small 
factor for a fixed total energy of the gas. Thus, we do not think that the results in our paper 
change significantly. Second, although recent works showed that jet matter is suggested to be 
electron-positron dominated at least at close to the core (e.g, Wardle et al. 1998), some authors 
proposed that protons constitute a part of the jet energy (e.g Mannheim 1998). If this is the case, 
the y value due to the jet matter is accordingly reduced, and the input energy epoch may not be 
constrained by SZ effect. 

Recently, Chandra observations found that there is no indication of shocks around several radio 
galaxies in the center of clusters (McNamara ct al. 2000; Fabian et al. 2000). These sources reside in 
such high pressure regions that the expansion speed of lobes is subsonic ( e.g. Fujita 2001). We have 
considered radio galaxies in proto-ICM which has not fully collapsed to have such high pressure, 
and then have shown that this assumption is not compatible with observational constraints of y 
and Einput- Thus Chandra findings do not alter our main conclusion. 

We thank useful comments of N. Sugiyama and M. Nagashima. We greatly appreciate B. Nath 
for reading this paper critically. 
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Fig. 1. — A contour map of the Compton y-parameter (dashed lines) and -Binput (sohd hnes). Each 
Une corresponds to y/10~^=l, 1.5, 3, 4, 5, 6 and £'input=0.14, 0.44, 0.74. 2.2, 3 keV, as indicated 
in the figure. A shaded region is allowed one by COBE observation (y < 1.5 x 10~^) and by the 
non-gravitational heating obtained by Lloyd-Davies et al. (2000) (0.44 ±0.3 keV). 
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Fig. 2. — Jet ignition epoch as a function of dark halo mass of a radio galaxy with various values of 
tgap- A shaded region corresponds to the allowed value obtained from the Figure 1 (tgap > 6.3 x 10^ 
years). Solid line (tgap = 0) indicates the typical formation epoch of objects of mass Mhaio- 



